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Abstract. Using numerical simulations we show how to realise an optical black hole 
laser, i.e. an amplifier formed by travelling refractive index perturbations arranged so 
as to trap light between a white and a black hole horizon. The simulations highlight 
the main features of these lasers: the growth inside the cavity of positive and negative 
frequency modes accompanied by a weaker emission of modes that occurs in periodic 
bursts corresponding to the cavity round trips of the trapped modes. We then highlight 
a new regime in which the trapped mode spectra broaden until the zero-frequency 
points on the dispersion curve are reached. Amplification at the horizon is highest 
for zero-frequencies, therefore leading to a strong modification of the structure of the 
trapped light. For sufficiently long propagation times, lasing ensues only at the zero- 
frequency modes. 
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1. Introduction 

Since Unruh's first suggestion that flowing media can be used as analogues for gravity pQ 
and in particular for studying certain physical phenomena usually associated to gravita- 
tional event horizons, e.g. Hawking radiation, the field has seen a steady increase in the 
number of proposed physical systems in which to observe such effects [21 13] . Originally, 
analogue gravity was proposed in systems that physically display a flowing medium 
that reproduces the flow of space close to a black hole: a flowing fluid that exhibits 
a gradient from subsonic to supersonic flow presents a point within the gradient such 
that the flow speed equals the speed of acoustic waves. This is a point of non-return 
for acoustic waves trying to propagate against the flow and is the analogue of an event 
horizon, causally disconnecting two separate regions within the medium. Such horizons 
are much more than mere toys or superficial analogies: they can truly reproduce the 
kinematics of gravitational event horizons. This implies that a direct analogue of the 
amplification of vacuum fluctuations and consequent emission of particles away from 
the horizon, known as Hawking radiation, should be visible. This idea has of course 
attracted significant attention along with the proposal of various systems that may en- 
hance the Hawking emission [3] . The main point is that Hawking emission is predicted 
to occur with a blackbody distribution with a temperature T connected to the horizon 
through the so-called surface gravity, k, i.e. the gradient of the gravitational field across 
the horizon [H [5] . In analogue systems, the surface gravity is determined by the gradient 
of the medium flow across the horizon and the combination of steep gradients along with 
very cold environments (in order to reduce the background thermal noise) is expected 
to lead to the observation of the spontaneous Hawking emission. Recently Philbin et 
al. proposed an optical analogue for horizons in which the steep gradients are ensured 
by shock front formation in intense optical pulses and the reduction of the background 
thermal noise is ensured by the fact that at optical frequencies, room temperature ther- 
mal noise is completely negligible [6J. 

In analogue models, black holes and their time-reversed realisation, i.e. white holes [7] 
play equivalently important roles. Whilst gravitational white holes are expected to be 
extremely rare, or even non-existent, analogue models abound with ways to generate 
such objects. Most notably, recent experiments performed in the presence of a white 
hole horizon in flowing water highlighted a stimulated frequency conversion process that 
can be considered as the classical limit of Hawking emission 0, E] whilst an optically in- 
duced white hole horizon was used to observe both the classical shifting of light frequency 
[13 QUI El H2] and the spontaneous emission of photons [131 EH] that has stimulated an 
ongoing discussion regarding the precise origin of the emission [151 [13 EE] ■ 
A remarkable prediction that emerged from analogue gravity studies was the predic- 
tion that a combination of a black hole and white hole event horizon would lead to a 
feed-back amplification process very similar to a laser where radiation is continuously 
bounced back and forth between the two horizons and hence exponentially amplified 

HB|. 
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Black hole lasers were first proposed proposed for systems that exhibit a "superluminal" 
dispersion relation (i.e. in which the group velocity increases with frequency) such as, for 
example, Bose-Einstein condensates (BECs). Black hole lasing has subsequently been 
studied both theoretically and numerically in Bose-Einstein-Condensate systems and 
similar settings [TSJ [IJJJ [201 121] • In this work we present a numerical study of black hole 
lasing in optical systems. Our numerics account for realistic parameters of the medium 
and highlight a laser amplification mechanism that is similar to that studied in BECs. 
However, we also show that as evolution ensues, the laser develops low frequency com- 
ponents that eventually seed amplification at what we call the "zero-frequency" modes. 
This amplification regime is exponentially favoured over the traditional black hole laser 
regime and leads to a significant increase of the output radiation. We then highlight 
and study how the low (close to zero) frequency region of the spectrum determines the 
properties of the high frequency emission, a feature that is expected to be generic to all 
analogue horizon systems. 

2. Optical horizon analogues. 

Optically induced horizon analogues are based on the use of nonlinear optics to create 
an effective moving medium [HI HH [22] • An intense laser pulse with carrier frequency 
Up, typically in the near-infrared or visible region, is focused into a dielectric medium 
and through the nonlinear Kerr effect it will excite a material polarisation response that 
contains two terms. One term oscillates at frequency 2u p and may excite photons from 
the vacuum state under the condition that energy and momentum conservation relations 
are satisfied. Unless particular care is taken in order to enforce such conditions, this 
so-called four-wave-mixing term (four elementary excitations are involved, two from the 
input pulse and two from the vacuum state) will be very strongly suppressed [23] • There 
is a second term that follows the envelope of the input pulse, i.e. it responds to the DC 
component of the electromagnetic field and is described in terms of a variation of the 
medium refractive index, n = n + n 2 I {(,)■, where n is the background refractive index 
(no ~ 1.45 in most glasses) and n 2 is the nonlinear refractive, also called "Kerr" index 
of the medium and ( = z — vt is the local lognitudinal cordinate. is the propagating 
intensity (envelope) profile of the laser pulse. In most common media n 2 ~ 10~ 16 W/cm 2 
so that the maximum index variation, 5n ma x ~ 0.005 — 0.0001, depending on the laser 
pulse intensity. With other media it may be possible to obtain higher 5n, e.g. some 
so-called soft glasses or media such as graphene are expected to give values larger by an 
order of magnitude or more. 

Light travels slower in high refractive index regions and faster in low refractive regions, 
so the moving laser pulse creates a moving 5n that locally slows light down. The effective 
space-time metric associated to the moving Sn is the so-called Gordon metric (2U [25] 
and may be recast in a form that tightly resembles the Painleve-Gullstrand metric for 
black holes [26] [23 [13] • In other words, in the reference frame comoving with the 5n, 
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Figure 1. Numerical simulations of the interaction of a laser pulse with a white 
hole horizon, (a) Amplitude profile, in logarithmic scale, of the laser pulse reflecting 
on the horizon (dashed line), (b) Evolution (along the propagation direction, z) of 
the spectrum, in logarithmic scale, during reflection from the horizon, (c) Schematic 
overview of the the geometry of a generic refractive index perturbation used in this 
work and location of the black and white hole horizons, (d) Dispersion curve used in 
this work with indicated the positions of the IN (input), P (positive), N (negative), 
Z (zero-frequency) and T (transmitted) modes. All modes are determined by the 
condition that frequency is conserved in the comoving frame, i.e. by the intersections 
of horizontal lines passing through the IN mode frequency, ±wjj^. 

space is flowing with velocity V = j 2 v(n 2 — l)/n: on the leading edge of the perturbation 
space flows inwards with increasing velocity until V — c (i.e. v = c/n) and a black hole 
horizon is formed. On the trailing edge, the time reversed situation is verified with 
space flowing outwards at a decreasing velocity and when V = c a white hole horizon is 
formed (see also [22J). 

In the presence of dispersion, i.e. of a frequency dependence of the medium background 
refractive index, one must take care in defining the exact nature of the horizon due 
to the fact that the phase and group velocities differ. In the following we will always 
consider the case in which the horizon is a blocking point for the group velocity of light, 
in line with the generally accepted opinion that this kind of horizon is the analogue for 
photons (or waves in a flowing medium) of a gravitational event horizon. 

3. Horizons in dispersive media. 

We now study how an input wave behaves when interacting with such horizons. 
This behaviour is analysed by performing numerical simulations based on Maxwell's 
equations, discretized and solved following the so-called Pseudo-Spectral-Space-Domain 
(PSSD) technique |28j that is a variation of the more common FDTD technique (29] 
and allows to include arbitrary material dispersion in a very simple manner. In more 
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detail, we consider 1-dimensional propagation modelled by the Maxwell equations for the 
transverse x and y components of the electric (E x ) and magnetic fields (H y ), respectively: 

d z E x = - n Q d t H y (1) 
d z H y = - d t D x (2) 

The electric displacement field is numerically evaluated as D x (t) = J r ^ 1 {e(u)J r [(l + 
e(z — vt))E x (t)]}, where J 7 and J 7 " 1 indicate the Fourier and inverse Fourier transforms, 
e(l>j) is the dispersive part of the medium refractive index, no = \fs. The moving 
refractive perturbation is described by a super-gaussian function e{z— vt) = Eo exp[— (z— 
vt) 2m /o~ 2m }, where m is the super-gaussian order and controls the perturbation steepness 
and a gives the perturbation width, typically fixed at 50 [im. The black hole cavity is 
constructed by using two such perturbations, displaced along the C-axis by a distance 
L that fixes the cavity length (see Fig. [3]). The input seed laser pulse is made to 
start directly inside the cavity and has the functional form E(z — vt) = exp[— (z — 
vt) 2 /a 2 E ] cos(kz — ut), where the initial pulse width <je is always taken to be shorter 
than the cavity length L, e.g. o~e ~ 3 — 6 urn and L ~ 20 — 80 /im. The photon numbers 
generated at each frequency are evaluated as S(u,z)/u, where S(cj,z) is the spectral 
intensity evaluated at each propagation distance, z. Total output photon numbers are 
obtained by integrating this quantity over frequency. The code was implemented with 
a dispersion relation shown by the solid line in Figure [IJd) in (u/, u) coordinates, where 
uj' = bj — v k with k = (u/c)n(u), is the frequency in the comoving reference frame. 
As can be seen, we are using a simplified dispersion curve that does not present any of 
the typical resonances present in most media such as glass. However, it was chosen to 
qualitatively reproduce the dispersion of diamond, a material that is currently becoming 
rather popular for a variety of reasons, including its remarkable transparency range due 
precisely to the absence of material resonances from sub-THz wavelengths up to the near- 
UV. Most importantly, such a choice allows us to study the generic physics underlying 
the black hole laser without additional complications due to specific structures within 
the dispersion relation. 

Figures [T](a) and (b) illustrate the behaviour at a white hole horizon: the input mode 
IN is chosen with a low frequency (3 «m wavelength) such that it travels with positive 
group velocity in the comoving frame (the IN mode identified on the dispersion curve 
in (d) has positive gradient) and therefore catches up with the travelling 5n. Upon 
reaching the horizon, light is blocked and transformed in to two new modes, one with 
positive frequency, P, and one with negative frequency, N. These both have negative 
group velocity in the comoving frame and are therefore reflected away from the horizon 
(see (a)). These two modes are the Hawking modes generated at the horizon. Most 
importantly we have explicitly verified that, in agreement with the interpretation in 
terms of Hawking emission, these two modes are such that the difference of their 
squared norms (i.e. their photon numbers) is equal to one, \P\ 2 — \N\ 2 = 1 and the 
ratio of the squared norms decays exponentially with increasing comoving frequency, 
liVlVl^l 2 oc exp(-Cw'), where C = l(2nc)/(j 2 v 2 )}l/(al5n/d() is the decay constant 
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Figure 2. Evolution of a laser pulse during the first few bounces from a white- 
black hole cavity (logarithmic scale), (a) Evolution along the propagation direction, 
z, of the amplitude in the comoving frame as the pulse bounces back and forth, (b) 
Evolution of the pulse spectrum, (c) Schematic representation of the cavity behaviour 
and indication of the various modes that are produced during each round-trip. 



and d5n/d( is the perturbation gradient evaluated at the horizon [T3j H%] . These two 
relations imply that the horizon emission in the P-mode exhibits a planckian black- 
body spectrum as a function of frequency. Most importantly for the black hole laser, 
the first relation also implies that the total photon number in reflection from the white 
hole horizon \P\ 2 + |iV| 2 > 1, i.e. the white hole horizon acts as an amplifying mirror. 
For small frequencies the photon gain \P\ 2 + |iV| 2 scales as 1/Cu', i.e. the steeper the 
gradient the larger the horizon gain will be. This simple recipe is important in order to 
understand and optimise black hole lasers. 

The behaviour of the black hole horizon can be inferred from that of the white hole as 
one is the time-inverted version of the other. This implies that a single mode impinging 
on the black hole horizon will not lead to negative mode generation or more in general 
to positive/negative mode mixing. The P mode will simply convert to the IN mode and 
reflect away from the black hole horizon. The N mode on the other hand will convert to 
a second red-shifted mode that still has negative group velocity and is thus transmitted 
through the horizon. This transmitted mode is indicated with "T" in Fig. |T|(d) (the 
dashed dispersion curve is relative to the region inside the black hole). 
Mode mixing and amplification at the black hole horizon will occur if both positive 
and negative modes, e.g. the same modes generated at the white hole horizon shown 
in Fig. pja), are sent on to the horizon as initial input conditions. In this case we 
indeed have the time reversed situation of the white hole horizon. However, as pointed 
out in [12], this amplification now relies critically on the phases between the P and N 
modes. For the case of large bandwidth pulses such as those used here, some frequency 
components will have the correct phase for amplification whilst others will have the 
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Figure 3. Possible configurations of the optical black hole laser in the presence of 
subluminal dispersion. The seed laser pulse is injected directly in the cavity and then 
propagates first towards the white hole horizon (as indicated by the arrows). 

wrong phase. Moreover, it is extremely difficult to purposely control these phases due 
to the varying and relatively strong dispersion (i.e. accumulated phases in traversing 
the cavity) across the pulse spectrum. We find that on average the P and N modes 
meeting at the black hole horizon lead to amplification, although this amplification is 
typically smaller than that occurring at the white hole horizon. 

In the black hole laser configuration discussed below, the white hole horizon will be 
approached by a single IN mode whilst the black hole horizon will always be approached 
with a simultaneous combination of both P and N modes, i.e. the modes generated at 
the white hole horizon. Therefore, both horizons will act as amplifying mirrors. 

4. The optical black hole laser 

On the basis of this reasoning we can construct an intuitive picture of the behaviour 
of a black hole laser. The black hole laser is formed by placing a black hole horizon 
and a white hole horizon in close vicinity such that waves are trapped between the two. 
In order for this to happen care must be taken in choosing the correct arrangement 
and this depends on the dispersion. In previous studies based on phonon oscillations in 
Bose-Einstein-Condensates, a so-called superluminal dispersion was considered, i.e. the 
dispersion curve had positive curvature [IBJ [HJJ [201 121] • In our case we have a dispersion 
curve with a negative curvature and so-called subluminal dispersion. This implies that 
in the laboratory reference frame, the group velocity of light is always smaller than the 
phase velocity and also that it decreases with increasing frequency. Therefore, based on 
the results in Fig. [TJ if we want to trap light between two horizons then the white hole 
horizon must be placed before (at smaller () the black hole horizon. 
We illustrate this situation in Fig. [2] that shows a numerical simulation of the first few 
cavity round-trips of a laser pulse trapped between two horizons. The spectrum evolves 
periodically at each reflection from the horizon mirrors. We clearly observe the same P 
and N modes seen in Fig. [T] at the white hole horizon. At the black hole horizon the 
spectrum returns to the original input frequency with the addition a second very weak, 
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Figure 4. Evolution of the laser spectrum (in log. scale) and electric field amplitude 
at the simulation output are shown for two different cases: (a)-(b) A£ = 0.3 /im, and 
(c)-(d) A£ = 3 /im. All other parameters are identical for both simulations: input 
wavelength 2 /im, cavity length L=20 /im, <5n m ax =0.1. 



less red-shifted mode around u — 1.5 x 10 15 rad/s. This is the "T" mode mentioned in 
the previous section [see also Fig. []Jd)] that originates from the partial transmission of 
the N mode through the black hole horizon and is ejected from the cavity. This mode 
is also visible in Fig. [2]^a) although due to its very low intensity it does not play an 
important role in the laser dynamics reported in this work. A schematic representation 
of all the laser modes is shown in Fig. |2](c) that summarises the behaviour of the light 
pulse within the cavity and the interaction with the two horizon mirrors. 
We note that it is possible to achieve a black hole laser in two ways as shown in Fig. [3j 

(a) we may use a negative 5n, similar to a potential well which will trap light. This 
is somewhat difficult to achieve of we assume that the 5n is generated by a laser pulse 
through the Kerr effect due to the fact that most media have positive ri2- Alternatively, 

(b) we may obtain the same effect by using two positive-valued perturbations generated 
by two independent laser pulses that are placed a certain distance apart. 

As an illustrative case, we show in Fig. [4] the numerical results for a black hole laser 
that is obtained using a negative 5n (5n m ax = —0.1) in which the black and white hole 
horizons are separated by 4.5 /im. The input mode is a Gaussian shaped pulse, with 
3 /im carrier wavelength and pulse length of 2 cycles. Figure |4]^a) shows the electric 
field in the vicinity of the black hole laser after a propagation distance of 650 fim. The 
vertical dashed lines indicated the position of the cavity that is travelling from left to 
right. The input mode, originally localised only within the cavity has led to lasing in the 
form of bursts of pulses tunnelling and exiting through the white hole horizon (not to 
be confused with the T-mode that would exit the black-hole horizon yet is too weak to 
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Figure 5. Evolution of the normalised laser spectrum (in logarithmic scale, 3 decades) 
and normalised electric field amplitude (linear scale, from -1 to +1) - profiles are 
shown for increasing propagation distances as indicated in the graphs. Each 0.1 mm 
of propagation corresponds very closely (in this simulation) to one cavity round-trip. 
Parameters used in the numerical simulations are: input wavelength 2 /im, A£ = 3 
jUm, cavity length L=80 fim, <5n m ax = 0.01. The shaded regions in (c) indicate the 
position and (rescaled) shape of the two 5n perturbations forming the cavity. 



be seen in this simulation). These bursts are composed of a sum of both P and N modes 
(hence the beat signal that modulates the individual bursts) and occur once every cavity 
round-trip, indicated with n=l,2,3,4 in the figure. These bursts are strongly increasing 
in amplitude at each round trip, a clear indication of the gain provided by the black hole 
laser. We note that here the modes leaking out of the cavity are doing so as a result of 
a tunnelling effect: the refractive index increase from the background to the maximum 
value (denoted by A( in Fig. |3[a)) at the white hole horizon occurs over a distance that 
is shorter than the optical wavelength. The modes therefore tunnel out of the hole - 
this tunnelling may be controlled by reducing the gradient of the refractive index profile 
and indeed, for refractive index changes that occur over scales larger than the mode 
wavelength, the outgoing leakage is suppressed and the modes are completely trapped 
within the cavity (see below). Figure [4](b) shows the spectral evolution corresponding to 
Fig. ^ a): the periodic back and forth conversion between in the IN and P/N modes at 
each cavity round trip can be clearly seen with a gain of more than 5 decades after only 
n=4 round trips. These findings are qualitatively similar to those reported for black hole 
lasers in Bose-Einstein condensates |19j . However, such a setting will be experimentally 
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Figure 6. (a) Evolution of the photon number N (in logarithmic scale) inside the black 
hole laser cavity for a Sn amplitude of 0.01 and two different configurations shown in 
Fig. [3] (solid line - negative Sn, dashed line - positive Sn). (b) Photon number insied 
the laser cavity for various maximum Sn amplitudes, indicated in the graph. 



challenging both in virtue of the large Sn (although similar results were obtained with 
smaller Sn, at the cost of correspondingly longer propagation distances) and of the 
extremely short distances at play, i.e. the short separation between the horizons and 
sharp A( increase at the horizon from the background index to the maximum Sn. 
In Figs. |4]^c) and (d) we show the same simulation with a transition A£ that has been 
reduced to experimentally reasonable values, i.e. A( ~ 3 /im. As can be seen, the high 
frequency modes no longer leak out of the white hole horizon. This allows for a higher 
gain within the cavity itself, although this is largely counterbalanced by the detrimental 
effect of reducing the Sn gradient as this reduces the amplification at each reflection. 
Nevertheless, in this situation a stronger build-up of negative frequencies occurs inside 
the cavity, which partly leak out from the black hole horizon at each round trip in the 
form of a more intense and now visible "T" mode, as discussed in Fig. [l](d). The T mode 
is emitted in bursts corresponding to the cavity round-trips (indicated with n=l,2,3,4 
in the figure). 

Finally, we note that in the very last stages of this simulation we can clearly see that 
the spectrum is starting to exhibit two distinct and new features: a strongly red-shifted 
component close to 00 = and a component, indicated with Z that lies in between the 
P and N modes. In the following, we focus attention on these new modes and highlight 
a new mechanism by which the pulse spectrum is broadened at each bounce until far- 
infrared wavelengths are generated. These then seed and sustain amplification at the 
zero-frequency points of the dispersion curve (indicated with filled squares in Fig. [ltd)), 




which have a gain favoured by the aforementioned 1/u/ dependence. This results in 
remarkably stronger gain and significant lasing even in conditions that could be realised 
in experiments. 

4-1. Zero-frequency mode amplification 

Figures |5^a) and (c) show results from a simulation in which the first several cavity 
round trips are shown. The horizon separation is 80 /zm, with 2 fim input wavelength, 
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A( = 3 /im, 5n = 0.01 - the curves are normalised and displaced horizontally for viewing 
purposes. The optical modes bounce back and forth but without any apparent overall 
increase in energy for nearly the first mm of propagation. At z=0.8 mm, a series of 
relatively weak bursts are clearly visible corresponding to T mode emission from the 
black hole horizon at each round trip. Figure [6^a) shows the total photon number in- 
side the cavity. During the first ~ 1 mm or propagation, the energy inside the cavity 
oscillates. These initial oscillations do not appear to be a numerical artefact and seem 
to correspond to real light pulse dynamics within the cavity. They most likely arise 
from a beating between the different modes that are concomitantly excited within the 
cavity and that therefore interfere with each other giving rise to amplitude oscillations, 
in a similar fashion to the build-up dynamics that may be observed in traditional laser 
cavities. After this short transitory period, a different regime appears in which the 
energy oscillations are washed out by a smooth and strong exponential increase of the 
total energy in the cavity. The origin of this change in behaviour can be understood 
from Fig. [5|a). During the first mm of evolution, light bounces back and forth but is 
gradually spread out to very low frequencies. This low frequency light has two effects: 
(i) due to its long wavelength it will start to completely fill the cavity and (ii) it will 
seed mode conversion to high frequency modes that lie close to or at the intersection of 
the dispersion curve with the horizontal axis, i.e. it will excite modes that have uo' — > 0, 
indicated with filled squares in Fig. [5tb) (we name the high frequency intersection as the 
"Z" mode). In order to understand why these zero- frequency modes appear we recall 
that, as a consequence of the planckian distribution of Hawking emission, amplification 
at the horizon diverges as 1/u/ for small u'. Therefore, the laser gain experienced by 
these very low frequency modes is remarkably stronger with respect to all other modes 
and they will therefore quickly take over: high-gain lasing at the "zero frequencies" will, 
and is expected to dominate the black hole laser scenario. 

In Fig. [6^b) we show the cavity photon number for identical operating conditions and 
with varying maximum 8n amplitude, as indicated in the figure. The gain (slope of the 
curve) is seen to decrease with decreasing 5n amplitude as a result of the fact that we 
have fixed the distance over which the perturbation switches on, A( = 3 fim. Therefore 
the horizon gradient (i.e. surface gravity) and hence also the horizon amplification de- 
creases with decreasing 5n max . Nevertheless, significant laser gain is observed even at 
low perturbation amplitudes that hold promise for future experiments. 
The question now arises as to how close to zero in the laboratory frame does the u — > 
zero-frequency mode have to be in order for this process to occur? The question is rele- 
vant as real experiments cannot be expected to rely on modes for example in the radio 
frequency domain. We therefore repeated the simulations by including a low frequency 
filter so as to completely suppress all radiation below a certain threshold value. Fig.[7]^a) 
shows the results for a 5n amplitude of 0.1 and filters as indicated in the caption. As 
may be seen, a filter placed at a wavelength of 20 /im hardly modified the overall gain 
and a reduction of ~50% is obtained only when the filter is placed at a wavelength of 8 
jim.. In Fig. [Wb) we repeat the simulations for an experimentally realistic 5n amplitude 




Figure 7. Evolution of the photon number inside the cavity for (a) Sn = —0.1 and 
(b) Sn — —0.005. The simulations were performed by including a spectral filter that 
cuts all frequencies below a certain value. The shortest propagating wavelengths are 
(a) thick solid line - no filter, thin line - 20 /im, dashed line - 8 /im, dotted line - 4 /im. 
In (b) thick solid line - no filter, thin line - 60 /im, dashed line - 10 um. 



of 0.005. A filter placed at 60 /im does not modify the laser gain and a 50% cut is 
observed with a filter placed at 10 /im. This is extremely promising. Indeed, many 
media are transparent in to the terahertz region (30-100 /im). For example diamond is 
transparent up to roughly 100 /im. This would therefore seem to indicate that this kind 
of novel amplification process could be observed in real settings. 



5. Conclusions 



We have numerically analysed the possibility of observing black hole lasing using laser- 
pulse induced (or optical) horizons. The general behaviour of the black-white hole 
horizon cavity is very similar to the behaviour of similar cavities studied theoretically 
and numerically in Bose-Einstein condensates. Our numerical simulations have been 
carried out over longer propagation distances (i.e. larger number of cavity round-trips) 
with respect to previous studies and highlight the onset of a regime in which the cavity 
leads to a broadening of the mode spectra until the two comoving zero-frequency modes 
are excited (in the laboratory frame, one mode lies at zero-frequency, the other lies 
in the visible or ultraviolet region, depending on the specific dispersion relation of the 
medium). Due to the 1/u/ dependence of the horizon amplification, the zero-frequency 
modes oscillate in the cavity with a much higher gain. These findings are summarised 
in the accompanying video animation that shows the evolution of the typical situation 
studied here: the initial optical pulse seeded in to the cavity is actually much shorter 
than the cavity itself (differently from typical cases studied so far in which the cavity 
length was of the same order of the input wavelength) and bounces back and forth until 
the zero-frequency modes take over and completely fills the cavity. 
Our numerical simulations treat the case of a coherently seeded laser. Naturally, one 
would expect that if the laser were to be seeded by vacuum fluctuations, then most 
certainly laser oscillation would occur at the zero-frequency modes as these have the 
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highest gain and lowest lasing threshold. 

A remarkable and unique feature of black hole lasers is the coupling between the very 
low frequency and very high frequency components of the electromagnetic spectrum. 
However, more complicated dispersion relations that include also resonances at certain 
wavelengths are expected to complicate this picture although a full model that includes 
losses due to material absorption would be required in order to correctly model such a 
situation. 

These simulations underline the fact that an experimental demonstration of black hole 
lasing, although not simple may not be as far-fetched as one may think. Refractive 
index perturbations of the order of 0.005 (or smaller) will give rise to significant gain 
over distances of the order of 1 cm. The dispersion relation used in this work is very 
close to that of diamond, a material that is now widely used and can be shaped in to 
optical waveguides. Diamond exhibits not only a very simple dispersion curve but (as a 
consequence of this) also a remarkably wide transparency range that could easily sustain 
the zero-frequency amplification regime reported here. The laser cavity, i.e. the black 
and white hole horizons could be generated by focusing for example two independent 
and ultrashort laser pulses, e.g. at 800 nm wavelength: each pulse would create a 5n 
through the nonlinear Kerr effect and the relative delay between the two pulses would 
allow to control the laser cavity length. The main challenge would be to find a condition 
(laser wavelengths, durations, energies) such that the rising and falling edges of the laser 
pulses that form the horizons are maintained shorter than ~ 3 fim for significantly long 
propagation distances, e.g. > 1 — 2 mm. Future studies will hopefully unravel other 
settings or combinations of materials and wavelengths that may indeed lead to the first 
experimental black hole laser. In the meantime, this remains a fascinating scenario in 
which to combine technologies and ideas developed in the area of photonics to the study 
of flowing media and horizon physics. 
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